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Abstract 

Considering the fluctuations of spectral functions, we prove that if chaotic systems fulfill the 
Bohigas-Gianonni-Schmit (BGS) conjecture, which relates their spectral statistics to that of ran- 
dom matrices, therefore by virtue of Gutzwiller trace formula, the instability of classical periodic 
orbits is constrained. In particular for two-dimensional chaotic systems, the Lyapunov exponent X p 
of each periodic orbit p should be bigger than a minimum value A m i n > 0.850738. This opens the 
possibility of new constraints for a system to be fully chaotic, or the failure of the BGS conjecture. 

PACS numbers: 05.45.Mt, 03.65. Sq 



1 



I. INTRODUCTION 



A difficult task in classica 
Except for some few maps 



mechanics is to prove that a dynamical system is fully chaotic. 



2J and billiards 



3|-[5| where analytical proofs exist, the numer- 



ical exploration of phase-space is the only practical way to decide this question. However, 
small regions of regular motion around stable periodic points can escape the numerical 
analysis, invalidating the method. 

One way to approach this question is to study the periodic orbits of the system. Lin- 
earization of equations of motion near a periodic orbit p is characterized by the monodromy 
matrix M p . In an autonomous Hamiltonian system with d degrees of freedom, the dimension 
of the monodromy matrix is 2d — 2, and the modulus of the determinant is one. Eigenvalues 
of this matrix are usually written in the form e p , where X p are the so called Lyapunov 
exponents, for 1 < i < 2d — 2. If the modulus of at least one of the eigenvalues is bigger than 
one (or equivalently the real part of some X p is positive), then the periodic orbit is unstable 
and the dynamics in the vicinity is hyperbolic. In a fully chaotic system all periodic orbits 
should be unstable, and the dynamics is fully hyperbolic jg]. 

In particular, for two dimensional Hamiltonian systems, the eigenvalues of the mon- 
odromy matrix of a periodic orbit p take the simple form {e Xp , ±e~ Ap }. Hence full hyper- 
bolicity implies that for every periodic orbit p, \ should have a nonzero real part. 

On the other hand, quantization of Hamiltonian systems have shown that imprints of 
classical dynamics are present in both eigenfunctions {?| and eigenvalues. These relationships 
are the main subject of study of Quantum Chaos. Regarding the eigenvalues, two important 
conjectures have been well established: 

(i) Berry- Tabor (BT) conjecture js]: for generic regular systems the eigenvalues of its 
quantization are completely uncorrelated and their statistics is Poissonian. 

(ii) Bohigas-Giannoni-Schmit (BGS) conjecture [9|: in fully chaotic systems the eigen- 
values are correlated (in particular there exists level repulsion) and they show the same 
statistical properties as those of random matrix ensembles. Particularly, for systems with 
time reversal symmetry, the fluctuations are those of the Gaussian Orthogonal Ensemble 
(GOE). For systems without time reversal symmetry the comparison is done with the Gaus- 
sian Unitary Ensemble (GUE) 

The BT conjecture was corroborated numerically for several cases and analytically for 
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some of them. Nevertheless some counterexamples were found, and the problem remains in 
the definition of a generic regular system. The BGS conjecture was also checked numerically 



in many theoretical systems, specially quantum billiards lOj, LL3J arid even experimentally 
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13]. We stress that the connection between complex dynamics and Random Matrix 



Theory (RMT) started in the 1950 's, where it was applied to the statistics of resonances in 
atomic nucl ei |14J| . The conjecture was also proved in some specific (but less realistic) systems 
like graphs 115] and billiards on negative curvature surfaces [5J. Although for the latter case 

n 

some counterexamples were found depending on the tiling properties of the billiards |16j . 

In the way to prove these conjectures is essential to use the connection between classical 
dynamics and eigenvalues that exist using the semiclassical approximation. For sufficiently 
high energies the density of states 

p(E) = J2$(E-E J ) , (1) 

j 

splits into a smooth plus an oscillatory part p = ~p + p- The smooth part admits an expansion 
in powers of energy called the Weyl series, which depends on global properties of the system 
such as volume, surface and topology of the classical phase space, not regarding the integrable 
or chaotic dynamics. The oscillatory part writes as 

p(E) = 2 M E ) cos [S P (E)/H + u p ] , (2) 
p 

where the sum is made over all classical periodic orbits p of the system that exist at energy 
E, S p is the corresponding action of the periodic orbit of period r p = dS p /dE, and v p is 
a phase called the Maslov index which depends on the number of turns of the stable and 
unstable varieties along the periodic orbit. 

The value of the amplitude A p depends on the stability of the periodic orbit. Berry 
and Tabor 8] found a general formula for stable periodic orbits in terms of the actions and 
frequencies of the classical tori that support the orbit. If the orbit is unstable, the amplitude 



is given by the Gutzwiller formula 17] 



A = 'p (3) 

P 2nh r pV /|det(M p -/)| ' 

in terms of the period t p , the monodromy matrix M p (I is the identity matrix), and the 
repetition number r p when the orbit p is non primitive (r p = 1 for primitive periodic orbits). 



The semiclassical expansion ([2]) can at least explain qualitatively some aspects of the 
above mentioned conjectures. In general, around a stable periodic orbit there exists a 
continuous family of periodic orbits, all with similar periods and actions. The contribution 
of a family produces a large fluctuation on the density of states if we compare it with the 
contribution of an unstable periodic orbit, which is always isolated. Therefore, in integrable 
systems, large fluctuations of the energy levels are expected, allowing for the characteristic 
clustering of a Poissonian distribution predicted by BT conjecture. On the other hand, 
in fully chaotic systems, where all the periodic orbits are unstable, the spectrum has less 
fluctuations and more rigidity, a characteristic shared with random-matrix eigenvalues. 

We show in this work an even deeper relation between the BGS conjecture and the 
classical dynamics. For every 2-D fully chaotic system: 

• if there is at least one periodic orbit whose Lyapunov exponent \ p verifies < A p < A m i n , 
then the fluctuations of the quantum spectrum do not obey random matrix statistics; 

• if the quantum spectrum has random matrix fluctuations, then all classical periodic 
orbits have Lyapunov exponents X p > A m i n . 

The two propositions can be demonstrated in the same way, but they have different 
implications. The former implies that the BGS conjecture is false, and more conditions are 
needed for a chaotic system to have a spectrum with random-matrix fluctuations. The latter 
has implications on classical dynamics: not only all Lyapunov exponents should be positive, 
moreover they should have a minimum value. 

In order to prove these propositions, a general spectral function depending on two param- 
eters is defined in Section [III and its fluctuations are estimated semiclassically by periodic 
orbit theory. In section II III the variance of fluctuations is calculated in terms of the spectral 
form factor, deducing the constraint to the monodromy matrix for each periodic orbit. We 
analize the implications of this result, particularly on the Lyapunov exponents of periodic 
orbits of two-dimensional systems, giving some conclusions in Section IIVI 

II. SPECTRAL FUNCTIONS 

Let {0 < Ei < E2 < ■ ■ ■ Ej ■ • ■ } be a discrete spectrum coming from the quantization of a 
bounded Hamiltonian system. We define a spectral function depending on two real positive 
parameters with energy units 
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7=1 ^ ' 



(4) 



where f(x) is a continuous square integrable function. We assume for simplicity that f(x) 
is an even function, with a Fourier transform 



f(u) = / f(x) cos(lox) dx . (5) 

J — oo 

We ask for f(x) to be a localized function, as a Gaussian, that decays exponentially for |se| ^> 
1. If f(x) = log(l + e x )/(l-|-e x )-|-log(i-|-e~ x )/(i + e _:r ), the spectral function F corresponds 
to the grand canonical entropy of a non-interacting Fermi gas with chemical potential p, 
temperature T (measured in units of Boltzmann constant), confined in a potential with a 



single-particle spectrum {Ej} [l8j 



By definition, the sum (j4]) can be written as the energy integral of / times the density of 
states (OQ) 

/oo / E\ f°° 

f (^jr- J p(E) dE = T J f(x) p(p-Tx)dx , (6) 

where the integration variable E has been changed to x = (// — E)/T. Because f(x) is 
localized around x — 0, it is only necessary to know the density of states in a window of 
order T around the energy /i. Moreover, if \i is high in the spectrum and T fi, we can 
introduce the semiclassical expansion (j2j) for p. With these conditions, in the energy window 
where the integrand ([6]) is significant, all classical quantities needed in ([2]) remain almost 
constant. The main dependence on energy comes from the phase of the cosine, where the 
action S p is divided by h, arriving to 



h + p 



(7) 



p{p — Tx) ps p(/i) + 2 A p (/j,) cos 

p 

Introducing this approximation into the integral can be performed analytically and the 



spectral function F also splits into a smooth plus an oscillatory part 



F(p,T) = I Tp(p) , (8) 
F(p, T) = 2 T £ A p f cos (| + ^ , (9) 
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where Iq = f(x)dx = /(O) is a dimensionless number. In the latter expression, all 
classical quantities A p , S p , t p and v v are computed at energy E = /i. Formula (jHJ) is the main 
result of this section. Its structure is similar to semiclassical formula (j2J), with \x playing the 
role of the energy. However, if f{uS) decays exponentially for to 3> 1, the prefactor Tf(Tr p /h) 
introduces a cut-off for long periodic orbits. Therefore, series ([9]) for F can be convergent if 
parameter T is big enough to compensate the exponential proliferation of periodic orbits in 
a chaotic system (which is related to its topological entropy). 



III. VARIANCE OF THE FLUCTUATIONS 

For a fixed value of T, rapid fluctuations of F are observed as /i is varied. These fluctua- 
tions have some distribution and typical size depending on several parameters. To character- 
ize the fluctuations around some value fi it is necessary to define some averaging procedure. 
For any spectral function G that depends on the energy /i, the average is defined as 

<G(ai)> = -r- / GfrW , (10) 

where A/i is the size of the window where the fluctuations are analyzed. From Eq. (J5J), 
every periodic orbit defines an energy scale 5{i p = h/r p where the phase of the cosine varies 
2tt. The largest of these scales is E c = h/r m i n , coming from the shortest periodic orbit of 
period r min . Then, if A/i ^> E c we know that (cos(S , p (/i)/^)) = for every periodic orbit. 
On the other hand, if A/i \i the classical quantities are almost constant inside the energy 
window (see for a detailed discussion on the involved energy scales). From now on we 
consider the fluctuations only at sufficiently high energies /i where the two conditions are 
fulfilled 



E c <A/i</i, (11) 

implying (F) = 0. Using Eq. ([9]), the variance (F 2 ) is written as the average of a double 
sum over periodic orbits. The product of cosines may be expressed as one half the sum of 
the cosine of the sum of the actions and the sum of their difference. The term with the sum 
of actions vanishes by the average, arriving to 
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F 2 ) = 2T 2 /j2A p A p/ f(r p T/h)f(r pf T/h) cos ) • (12) 

\ p,p' I 
We can relate the variance with the semiclassical expression of the form factor (the Fourier 

transform of the two point correlation function of the density of states) [20j 

K(E, r) = h* (j2 A P A P , cos (^^) 6 (r - ^±^) ^ . (13) 
Given this expression, the variance of the fluctuations is written in the integral form 

( F2 )-2^l /VWMr. (14) 
All the dynamical information of the system is now contained in the form factor. It is 
also known that for times r < th, where th = h~p(E) is the Heisenberg time, the only 
contributions to the form factor come from terms where p' = p, and p' = p T if the orbit p 
has a time reversal partner p T . This is the so called diagonal approximation, and the form 
factor simplifies to 

K D (E, r)^h 2 J2 9 P A 2 p 8{t - r p ) . (15) 
p 

In systems with time-reversal symmetry g p = 2 for all periodic orbits (excepting those few 
ones that are self- retracing, where g p — 1). If time- reversal symmetry is broken, q„ = 1 for 



all periodic orbits. For chaotic systems, by Hannay-Ozorio de Almeida sum rule [2l|, this 
simplifies to K D (fi,r) = g r for r > r min (g = 2 for time-reversal systems, g — 1 if the 
symmetry is broken) . This approximation is valid for the full form factor up to a time much 
smaller than Heisenberg time. Thereafter, off-diagonal corrections become important and 
expression ( fl3l) converge to the random matrix form factor 



K goe (E,t) = < 



' 2r-rlog(2^ + l) , r < r H 
^T H -Tlog(g±m) , T>T H 



(16) 



r , r < r H 



K gve (E,t) = < 



(17) 



m , r > r H 
7 



GOE (GUE) should be used for chaotic systems with (without) time reversal symmetry. For 
t ^> th the form factor goes to a constant value th for both symmetries. 

If we consider the contribution to the variance of F given by all long periodic orbits 
t p > Ti, such that r min < r! < t h , we have 

rp2 POO 

V ( T ^ = 2^Wj f 2 (rT/h)K^,r)dr . (18) 

For this time regime the form factor should correspond to the random matrix prediction. 
Moreover, we know that K KMT (n, r) < Ko(fi, r) = gr, arriving to the inequality 

nT 2 r°° - 

V(n) < JL— / rf(rT/h)dr . (19) 



2ir 2 h 2 

Changing the integrating variable to u — tT/H the integral simplifies to 



< 7T1 / "f 2 (") d " » (2°) 
where U\ = r\Tjh. Because the integrand is always finite and positive, we can extend the 
lower limit up to zero 

/»oo 

V[Tl)< hl "/»d" = ^. (21) 
The integral converges to I±, assuming that f(u) decays rapidly enough, which imposes fur- 
ther restrictions on the possible primitive functions f(x) that can be used. The contribution 
to the variance of long periodic orbits has an upper limit that depends only on the Fourier 
transform of the function f(x). 

On the other hand, using Gutzwiller formula ([3]) for A p in chaotic systems, we can write 
the contribution to the variance from long periodic orbits in the semiclassical form 

Wn )-—Y s 9p ^ f 2(UJp) (22) 
V{Tl) -2n 2 ^ Ti r 2 \det(M p -I)\ ' (22) 

where u p = r p T/h. In this sum, all terms are positive. For some u> = u>* the function 
uj 2 f 2 (uj) has a finite maximum value I*. Therefore, if we consider the contribution of one 
periodic orbit p to the sum (I22j) . and we fix the parameter T to T* = u*h/r p , we show that 

,/( -) > 2^||^W-7)| ' (23) 
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Comparing this to the inequality (|2T|) . that holds for every T, we arrive to the condition 



9 ' <±h. (24) 



27r 2 r2|det(M p -/)| 2tt 2 " 
For the general case of primitive periodic orbits (r p = 1), and non self-retracing orbits 
(dp = 9)1 condition 021]) simplifies to 

|det(M p -/)| > a = — . (25) 

h 

This is the key result of the work, showing a constraint on the instability of periodic orbits 
in chaotic systems. For two-dimensional systems the Lyapunov exponents of an unstable 
periodic orbit are {X p , —X p } by the condition that | detM p | = 1. Therefore 

|det(M p - I)\ = 4sinh 2 (A p /2) , (26) 
and inequality ( 125]) can be inverted giving 

X p > 2 arcsinh(Va/2) . (27) 

This condition should be fulfill for every even symmetric function f(x) with finite Iq integral 
and finite I\ integral related to its Fourier transform. We face a difficult variational problem: 
to find the function that maximizes the value of a, in order to find the threshold value A m i n 
for the Lyapunov exponent in 2D systems. A first test with some elementary functions has 
given a maximum value a = 2/e = 0.73576 for the Gaussian function f(x) = e~ x2 . From 
this result we explored the family function f(x) = e - ^'", varying the exponent n. By a 
numerical analysis we found the maximum value a = 0.76847 for exponent n = 3, and using 
(]27p the highest limit to the Lyapunov exponent is at least 

X p > A min > 0.850738 .... (28) 

Some other families of functions were tested although none of them had given a higher value 
of a. 
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IV. CONCLUSION 



We have shown that if BGS conjecture is true, i.e. spectral fluctuations of fully chaotic 
systems follow RMT, then instability of classical periodic orbits has a constraint. Particu- 
larly in 2D systems Lyapunov exponents of long periodic orbits have a minimum value. 

As it was mentioned in Section m a physical example of a spectral function F(/j,, T) is 
the grand canonical entropy of a non-interacting fermion gas in a confining potential. For 



this detailed analysis of the behavior of ( F 2 j as a function of parameter T (the 

temperature) was done in [19] for gases with chaotic single-particle motion. For T<A the 
variance of entropy fluctuations grows linearly with T. Increasing T up to few mean level 
spacings, the variance saturates to the maximum value V ma3L = gl\j {2tt 2 ) given by RMT, 
see Eq. fl2T]) . When T is of order E c , the shortest periodic orbits become relevant and an 
exponential decreasing of the variance is expected. 

This analysis will be qualitatively similar for other localized functions f(x). Therefore, 
even if the proof of (125]) was made for long periodic orbits r p ^> r min , where usually the 
Lyapunov exponents are large, we conjecture that this condition can be extended to the 
shortest periodic orbits. Violation of (I2"5"|) for a short periodic orbit p would imply that the 
variance of the spectral function F(fi, T), after reaching the maximum value V max , instead of 
decreasing exponentially to zero, would continue to grow beyond the value given by Eq. ( 123]) 
at some value T* = u*h/r p . Numerical calculations in 19) for a Sinai-type billiard, and in 
22j for a fictitious system given by the non-trivial zeros of the Riemann Zeta function, have 
shown that the variance of entropy fluctuations, as a function of temperature, never become 
bigger than the limit imposed by RMT. We induce from these observations that monodromy 
matrix of all periodic orbits in these systems fulfill condition ( 125]) . 

Nevertheless, an extensive numerical analysis of other chaotic systems is needed to give 
more evidence to the conjecture for short periodic orbits. Or eventually a counterexample 
would keep the constrain |det(M p — I)\ > 0.76847 only for long periodic orbits, as it was 
proved in this work. 
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